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ABSTRACT 

We examine the stellar mass assembly in galaxy cluster cores using data from the 
Cluster Lensing and Supernova survey with Hubble (CLASH). We measure the growth 
of brightest cluster galaxy (BCG) stellar mass, the fraction of the total cluster light 
which is in the intracluster light (ICL) and the numbers of mergers that occur in the 
BCG over the redshift range of the sample, 0.18< 2 < 0.90. We find that BCGs grow 
in stellar mass by a factor of 1.4 on average from accretion of their companions, and 
this growth is reduced to a factor of 1.2 assuming 50% of the accreted stellar mass 
becomes ICL, in line with the predictions of simulations. We find that the ICL shows 
significant growth over this same redshift range, growing by a factor of of 4-5 in its 
contribution to the total cluster light. This result is in line with our previous findings 
for ICL at higher redshifts, however our measured growth is somewhat steeper than is 
predicted by simulations of ICL assembly. We find high mass companions and hence 
major merging (mergers with objects of masses ^1/2 of the BCG) to be very rare for 
our sample. We conclude that minor mergers (mergers with objects with masses < 
1/2 of the BCG) are the dominant process for stellar mass assembly at low redshifts, 
with the majority of the stellar mass from interactions ending up contributing to 
the ICL rather than building up the BCG. From a rough estimate of the stellar mass 
growth of the ICL we also conclude that the majority of the ICL stars must come from 
galaxies which fall from outside of the core of the cluster, as predicted by simulations. 
It appears that the growth of the ICL is the major evolution event in galaxy cluster 
cores during the second half of the lifetime of the Universe. 

Key words: galaxies: clusters: general - galaxies: clusters: intracluster medium - 
galaxies: interactions - galaxies: evolution - galaxies: elliptical and lenticular, cD 


1 INTRODUCTION 

Having formed from the largest density perturbations re¬ 
maining after the Big Bang, galaxy clusters give us insight 
into the conditions in the early Universe and can provide 
constraints for cosmology. 

The most massive galaxies in the Universe reside in the 
cores of galaxy clusters. Brightest cluster galaxies (BCGs) 
generally sit close to or at the centre of the gravitational 
potential of their host cluster. BCGs are thought to have 
formed and evolved by a hierarchical merging process, which 
would cause them to grow gradually over time. This grad¬ 
ual assembly is predicted by cosmological simulations of 


structure formation (e.g. |Dc Lucia fc Blaizotj 2007), how¬ 
ever observations of the stellar masses of BCGs show that 
they have assembled the majority of their mass fairly early 
on, with some studies finding as much as 95% of present 
day BCG masses being already in place by redshift 2 ~ 1 


(equivalent to about half-way back to the Big Bang; Collins 
|et al.|2009| |Stott et al.|2010' I. Subsequently a range of mass 
growth estimates have been reported over the timescale be¬ 
tween 0 < 2 < 1, with some studies showing more growth 
over this redshift range (e.g. Lidman et al.|[2012 Lin et al. 


2013) necessitating a less extreme and rapid early assembly. 


In any case the observed masses of BCGs at high redshifts 
are generally larger than are predicted by simulations. This 
apparent rapid and early assembly of massive galaxies in 
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cluster cores is a major challenge for models of galaxy for¬ 
mation and remains a topic of much debate. 

Galaxy cluster cores are observed to be pervaded by a 
diffuse, low brightness, stellar component known as intra¬ 
cluster light (ICL). The ICL is made up of stars which are 
thought to be gravitationally bound to the cluster potential 
rather than any specific galaxy. The origin of the ICL is still 
unclear at the present time, although it is generally thought 
that it can be produced by interactions between galaxies 
in clusters, such as tidal stripping and merging. In nearby 
clusters the ICL, whilst spread between the majority of the 
cluster galaxies, is often found to be centrally concentrated 
around the BCG (e.g. |Mihos et al.|2005||Rudick et al.|2011[ ). 
As such it seems that the evolution and assembly histories of 
BCGs and the ICL are probably intertwined. Observations 
and simulations of the ICL show a large growth between 


0 < z < 1 (e.g. Kric.k & Bernstein|2007 

Murante et al. 

2007j 

Rudic.k et al. 2011 Burke et al. |2012 

Contini et al. 

2014} 

and when combined with the expecte 
ers in cluster cores and onto BCGs (e 

1 numbers of merg- 

.g. Burke & Collins 

2013( Edwards & Patton|2012||Laporte et al.|2013; Lidman 

et al. 112013 ( this presents a scenario of the ICL being built 


up by interactions of galaxies with the BCG. This mecha¬ 
nism for ICL growth would allow BCG masses to remain 
relatively unchanged, or to only grow by a small amount 
over this redshift range if the majority of the stellar mass 
from mergers ends up in the ICL rather than centrally on 
the BCG. Simulations of ICL buildup and mergers in cluster 
cores show that a large fraction of merging stellar mass (be¬ 
tween 30-80%, see e.g. Conroy et al.|[2007l |Puchwein et al.| 


|2010| |Laporte et al.|2013 ) will end up residing in the diffuse 

halo of the BCG or ICL, however this scenario has not yet 
been observationally confirmed. 

In this paper we aim to examine the growth and buildup 
of BCGs and the ICL, and test the scenario of ICL growth 
through interactions between galaxies in the cluster core. 
In order to do this we use the X-ray selected Cluster Lens- 


ing and Supernova survey with Hubble (CLASH; Postman 


et al. 2012) cluster sample which covers the redshift range 


0.18 < z < 0.89. We measure the BCG stellar masses, the 
contributions of the ICL and BCG to the total cluster star 
light, and the number of mergers that the BCG may have ex¬ 
perienced over this redshift range. Our approach represents 
a self-consistent analysis of the mass assembly in galaxy 
clusters including BCG assembly, ICL and BCG compan¬ 
ion galaxies. 

The structure of this paper is as follows: Section [2] con¬ 
tains the details of the sample studied; in Section [3] the 
methods and corrections used are described; in Section [f] 
we present the results of our measurements; the results are 
discussed in Section [5] and conclusions drawn in Section [6j 
Throughout this paper we adopt a ACDM cosmology with 
Ho = 70 km s -1 Mpc -1 , Qm = 0.3, Ha = 0.7. All magnitudes 
quoted are in the AB system. 


2 DATA 


CLASH (Postman et al. 20121 contains 25 galaxy clusters 
in the redshift range 0.18 < z < 0.89 which were selected 
to have temperatures Tx >5 keV. Some 20/25 clusters were 
selected for their X-ray surface brightness symmetry indi¬ 


cating relaxed clusters, and the remaining 5 were specifically 
selected to have large Einstein radii. The clusters were se¬ 
lected with their lensing properties in mind and at least 18 
of the clusters show strong lensing arcs, however the X-ray 
symmetry selection criteria goes some way to ensure that 
these clusters are not preferentially aligned along the line- 
of-sight. We use 23 of these clusters which are described 
in Table |T| The clusters we exclude are MACS1931.8 and 
MACS0717.5 which both have bright stars near the cluster 
core, the light from which contaminates the ICL and overall 
cluster light in a way which we cannot accurately quantify 
or remove. CLASH contains 16 band HST photometry for 
each cluster, the data used here were downloaded in reduced 
form from the CLASH websittQ 

Since the ICL is extended and has low surface brightness 
it could be easily mistaken for sky background during data 
reduction and subtracted. Presotto et al. (2014} examine one 
of the CLASH clusters (MACSJ1206, z = 0.44) using multi¬ 
band SUBARU data, and as part of their data reduction 
they test that sky subtraction does not cause the ICL to 
be removed by comparing their Subaru reduction with the 
HST data used here. They find the ICL to be detected to the 
same surface brightness levels in both data sets. Since both 
data use independent sky subtraction methods this shows 
that the over-subtraction of ICL which is mistaken for sky 
is not an issue here. 


3 ANALYSIS 


Before we analyse the data we first remove any point sources 
and foreground or background galaxies from the images. 
Point sources and non-cluster member galaxies are identified 
using the object catalogues for the CLASH sample described 
in Jouvel et al. (2014}. Point sources are defined from the 
CLASH catalogues as objects having a “stellarity” value of 


0.1 or greater. In line with Jouvel et al. (20141, a stellarity 


less than 0.08 selects 86% of objects in any given CLASH 
cluster image - a value of 0.1 effectively identifies all the 
galaxies in an image. The point sources in the data were 
thus identified and were masked out using a circle of area 
1.5 times the area defined for that object in the catalogue. 
The data were then visually examined and any remaining 
obvious point sources were removed by hand. Cluster mem¬ 
ber galaxies are identified as any object whose photometric 


redshift from the photo-z catalogue (Jouvel et al. 2014) is 


within 3cr of the cluster redshift, where the 1 <t accuracy of 


the photometric redshifts is defined in Jouvel et al. (2014) 


as 0.04 x (1Y Zobject)• Non-cluster galaxies were also masked 
out using a circle of 1.5 times the area given in the cata¬ 
logue, in the case of gravitational lens arcs, the arcs were 
masked out by hand. 


3.1 SED fitting and stellar mass measurements 

To measure the stellar masses of the BCGs and galaxies 
within 50kpc of the BCG in the CLASH sample we fit 
their spectral energy distributions (SEDs) using photometry 


1 http: 

archive .stsci.edu/prepds/clash/ 
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Table 1. Clusters from CLASH used in this study. 


Cluster 


Redshift RA Dec Cluster mass M 200 

xlO ls M 0 


Cluster mass reference 
(from lensing measurements) 


Abell 383 

0.187 

02 

48 

03.36 

-03 

31 

44.7 

1.04+ 

0.07 

Abell 209 

0.206 

01 

31 

52.57 

-13 

36 

38.8 

1.17+ 

0.07 

Abell 1423 

0.213 

11 

57 

17.26 

33 

36 

37.4 

1.20+ 

0.59 

Abell 2261 

0.224 

17 

22 

27.25 

32 

07 

58.6 

1.76+ 

0.18 

RXJ2129+0005 

0.234 

21 

29 

39.94 

00 

05 

18.8 

0.73+ 

0.18 

Abell 611 

0.288 

08 

00 

56.83 

36 

03 

24.1 

1.03+ 

0.07 

MS 2137-2353 

0.313 

21 

40 

15.18 

-23 

39 

40.7 

1.26+ 

0.06 

RXJ1532+30 

0.345 

15 

32 

53.78 

30 

20 

58.7 

0.64+ 

0.09 

RXJ2248-4431 

0.348 

22 

48 

44.29 

-44 

31 

48.4 

1.40+ 

0.12 

MACSJ1115+01 

0.352 

11 

15 

52.05 

01 

29 

56.6 

1.13+ 

0.10 

MACSJ1720+35 

0.391 

17 

20 

16.95 

35 

36 

23.6 

0.88+ 

0.08 

MACSJ0416-24 

0.396 

04 

16 

09.39 

-24 

04 

03.9 

2.50+ 

0.50 

MACSJ0429-02 

0.399 

04 

29 

36.10 

-02 

53 

08.0 

0.96+ 

0.14 

MACSJ1206-08 

0.440 

12 

06 

12.28 

-08 

48 

02.4 

1.00+ 

0.11 

MACSJ0329-02 

0.450 

11 

57 

17.26 

-02 

11 

47.7 

0.86+ 

0.11 

RXJ1347-1145 

0.451 

13 

47 

30.59 

-11 

45 

10.1 

1.35+ 

0.19 

MACSJ1311-03 

0.494 

13 

11 

01.67 

-03 

10 

39.5 

0.53+ 

0.04 

MACS J1149+22 

0.544 

11 

49 

35.65 

22 

23 

55.0 

5.10+ 

1.90 

MACSJ1423+24 

0.545 

14 

23 

47.76 

24 

04 

40.5 

0.65+ 

0.11 

MACSJ2129-07 

0.570 

21 

29 

26.06 

-07 

41 

28.8 

3.50+ 

3.10 

MACSJ0647+70 

0.584 

06 

47 

50.03 

70 

14 

49.7 

6.80+ 

1.40 

MACSJ0744+39 

0.686 

07 

44 

52.80 

39 

27 

24.4 

0.79+ 

0.04 

CLJ1226+3332 

0.890 

12 

26 

58.37 

33 

32 

47.4 

1.72+ 

0.11 

CLASH ACS (for this study we 

used the F606W, 3.2 ICL Meas 
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F625W, F775W, F814W, and F850LP bands) and WFC3 
(F105W, F110W, F125W, F140W, F160W) images. Mag¬ 
nitudes were measured in each band using the MAG_AUT0 
method in SExtractor (Bcrtin & Arnouts 1996; similar to 
the method described in Stott et al.|2010'l, using consistent 


To measure the ICL we follow the method described in IBurkel 


et al. (2012) which is summarised here. The total cluster 


aperture sizes between wavebands for each galaxy, set by the 
F140W band. 


We fit the resulting SEDs using a similar methodology 
to that described in |Shapley et al. (|2005[ see also |Hilton| 


et al. 2010 Menanteau et al. 2012). We construct a grid 


light is measured by summing all the flux in the images af¬ 
ter non-cluster objects have been masked out in the manner 
described above. The fraction of the total cluster light con¬ 
tained in the ICL is simply estimated by summing all the 
light below a given surface brightness threshold and com¬ 
paring this to the total cluster light. In this paper we use a 
threshold of 25 mag/arcsec 2 in the rest-frame B-band which 
is comparable to that used in previous studies (e.g., |Feld- 


of Bruzual & Chariot (2003I solar metallicity models with 


exponentially declining star formation histories with 20 val¬ 
ues of t in the range 0.1 — 20 Gyr, and 172 ages in the 


range 0.001 — 13.0 Gyr. We adopt a Salpeter (19551 initial 


mass function (IMF), and model the effect of dust extinc¬ 
tion using the Calzetti et al. (2000) law, allowing E(B — V) 
to vary in the range 0.0 — 0.48 in steps of 0.04. We use x 2 
minimisation to find the best fit and determine the stellar- 
mass of each galaxy from the normalisation of the best fit 
model template to the observed SED. We restrict the range 
of model templates to those with ages less than the Uni¬ 
verse at the redshift of each cluster. Despite the fact that 
degeneracies exist between model parameters (e.g., age, r, 
E(B — V); see the discussion in Shapley et al. 2005), many 
authors have shown that stellar masses measured from SED 
fitting are robust to variations in model parameters (e.g. 
|Forster Schreiber et al.|2004| |Shapley et al.|2005| ) . We note 
that uncertainties in the IMF and the contribution of ther¬ 


Imeier et al.|2004||Krick et al.|2006|| Krick fe Bernstein|2007| 

|Burke et al'|2012||Presotto et al.|2014 |. Whilst this method 

is very simplistic it avoids the necessity of disentangling the 
low brightness light at the edges of galaxies from the faint 
ICL, which is impossible to do definitively without dynami¬ 
cal information about the stars in the cluster. This method 
also avoids the necessity of fitting both the ICL and BCG 


together with a surface brightness profile (e.g., see Zibetti 
etaL|2005l |Gonzalez et~aI1|2007| pUTa] |Toledo et al.||2011| ). 
The simplicity of this method also makes it very effective 
for comparing ICL results across different studies and data 


samples (e.g., above references). See Burke et al. (2012) for¬ 


mally pulsating AGB stars to infrared fluxes (e.g. Maraston 
|2005| |Conroy et al.|2009| ) is likely to lead to the stellar mass 
estimates only being accurate to within a factor of ~ 2. 


a more detailed discussion of this method. 

Due to the clusters being at a variety of redshifts, the 
rest frame B-band needs to be converted to the appropriate 
observed frame waveband from which the ICL is to be mea¬ 
sured. This is done using the standard redshift-wavelength 
equation A 0 f, s = A em it/(1 + z), the observed wavebands 
for each cluster are listed in Table [2] The rest-frame 25 
mag/arcsec 2 surface brightness limit also needs to be con¬ 
verted to an observed surface brightness limit for each in¬ 
dividual cluster by applying the redshift surface-brightness 
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dimming correction (2.51og(l + z) 3 ) and the K-correction 
for each cluster’s respective redshift. The K-correction cal¬ 
culated for these observations is a differential K-correction, 
dK(z), taking into account the difference in central wave¬ 
length of the rest-frame B-band and the de-redshifted ob¬ 
served frame wave-band. The dK values were calculated us¬ 
ing the stellar population synthesis models of |Bruzual fe| 
Chariot (2003) assuming an old stellar population with a 


formation redshift of Zf = 3 and a solar metallicity. This is 
obtained using the magnitudes derived from the |Bruzual fe| 
Chariot (20031 models as, 


dK(z ) = M ev (z) — M r f(z) + 2.5 log(l + z), 


(1) 


where the M ev term is from the observed HST filter, and the 
M r f term is from the rest frame B-band. The values for dK 
are shown in Table [2] The observed band surface brightness 
limit is thus given by, 

UobsC) = Urest(z) + 2.51og(l + z) 3 + dK(z), (2) 

(see H ogg et al.|2 002 |. The corrected values for the observed 
surface brightness threshold are given in Table [2] Exam¬ 
ples of the cluster cores with the ICL thresholds applied are 
shown in Figure [l] The 2cr surface brightness limit for the 
data is 26.3, 26.4, 25.7 and 25.5 mag/arcsec 2 for the F606W, 
F625W, F775W and F850LP filters respectively. Since the 
observed surface brightness thresholds for the F775W and 
F850LP filters are fainter than the sky limit we cannot mea¬ 
sure the ICL below 25 mag/arcsec 2 for the higher redshift 
clusters in CLASH. All the clusters observed in the F625W 
filter above a redshift of 2 = 0.4 are also excluded from the 
ICL analysis for this reason (see Table [2|. 

The error on our measurement of the ICL fraction of 
the total cluster light is estimated from the variance of the 
flux (f) in the ICL as, 




flCL-CTtotal \ 2 


( ftotal ) 2 

where a denotes the standard deviation. 


(3) 


3.3 Mergers onto the BCG 


In order to estimate the stellar mass growth of the BCG and 
ICL we estimate the number of mergers which will occur 
with the BCG. To do this we follow the method described 
in |Burke fe Colliris] ( |2013| ) (hereafter BC13) of counting the 
number of companions to the BCG (N c ) within a fixed aper¬ 
ture. For comparison with previous studies we use an aper¬ 
ture of 50 kpc, centred on the BCG. The BCG companions 
were detected using SExtractor ( Bertin fe Arnouts|1996 1. 

We carried out completeness tests for the detection of 
BCG companions by inserting fake model galaxies, with (del 


Vaucouleurs (19481 profiles and half-light radii of 5 kpc (a 
similar size to the smaller BCG companions in the data), 
into the cluster images at random positions within a 100 kpc 
radius of the BCG. We then ran SExtractor over the im¬ 
ages and counted the number of model galaxies that were 
detected. We decreased the brightness of the model galaxies 
inserted until the detection completeness fell below 90%. We 
found our detection of companions to be complete to a min¬ 
imum brightnesses of 1:1000 of the BCG. All companions 
subsequently counted in our results have brightnesses above 


Table 3. Dynamical friction timescales for typical masses of com¬ 
panions measured from SED fitting assuming a mass-to-light ratio 
of M/L=4, ry=550 krn s’ 1 and an initial distance from the BCG 
of r=50 kpc. 


Stellar mass ratio 

Tfric (Gyr) 

1:2 

0.38 

1:5 

0.96 

1:10 

1.93 

1:20 

3.86 

1:50 

9.65 

1:100 

19.29 

1:1000 

192.93 


this limit; any companions fainter than this limit are ex¬ 
cluded from our analysis. The companions counted by SEx¬ 
tractor were also checked by eye to ensure that bright com¬ 
panions near to the BCG were not accidentally excluded and 
close pairs of objects were not miscounted. 

In order to be able to say whether or not a BCG com¬ 
panion will merge with the BCG in the time between the 
cluster redshift and the present we need to know its merging 
timescale. As described in BC13, we estimate the amount of 
time taken for a companion on a circular orbit to be accreted 
onto the BCG using the dynamical friction timescale, T/hc 
(G yrs), which is given by, 


Tfric = 3.81 x 10' 


,5 Cvc 

M c ' 


(4) 


Here r is the separation of the galaxies in kpc, v c is the 
circular velocity of the galaxy cluster in kms -1 , and Me is 
the total mass of the companion in solar masses (for a more 
thorough description see |Binney fe Tremaine||l987| and for 
a more detailed description of its use in this instance see 
BC13). 

We estimate the stellar masses of the BCGs and their 
companions by fitting their SED in the same manner as 
described above. In order to calculate their total masses 
(M c ) we then assume a mass-to-light ratio of M/L=4 (in 
the rest-frame B-band), in line with BC13 and references 
therein. In the absence of dynamical information, the cir¬ 
cular velocity of the companions is also estimated in the 
same way as BC13. Typical values of circular velocities in 
galaxy cluster cores range between v c ~ 350 kins -1 and 
v c ~ 700 kms - (see BC13). In order to calculate Tf we 
use an average of these and estimate the circular velocity 
to be v c = 550 kms - . These values for M/L and v c com¬ 
bined with the stellar masses measured from the SED fitting 
give the typical dynamical friction timescales, Tf r i c , and are 
summarised in Table [3] As is clear in Table [3] the least mas¬ 
sive BCG companions have merging timescales in excess of 
the Hubble time. The dynamical fraction timescale was cal¬ 
culated for each BCG companion individually, whose mag¬ 
nitude is above the completeness limit. Only companions 
which have dynamical friction timescales less than the time 
between the observed cluster redshift and the present day 
are counted as contributing to the BCG mass growth in the 
subsequent analysis. 

The errors on the numbers of BCG companions are cal¬ 
culated using the photometric redshift accuracy quoted in 


Jouvel et al. (2014). Jouvel et al. hnd 96% of the CLASH 
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Table 2. The corrections applied for ICL measurement and the observed surface brightness threshold equivalent to the rest frame limit 
of fiB =25 mag/arcsec 2 . 


Cluster 

z Observed waveband 

for ICL measurement 

M ev (z) — M r f(z) dK Surface brightness 

dimming correction 
2.5 log(l + z) 3 

Corrected surface brightness 
equivalent to 25 mag/arcsec 2 
at z = z (cluster) 

Abell 383 

0.187 

F606W 

-0.68 

-0.49 

0.56 

25.06 

Abell 209 

0.206 

F606W 

-0.62 

-0.41 

0.61 

25.20 

Abell 1423 

0.213 

F606W 

-0.55 

-0.34 

0.63 

25.29 

Abell 2261 

0.224 

F606W 

-0.55 

-0.33 

0.66 

25.33 

RXJ2129+0005 

0.234 

F606W 

-0.44 

-0.21 

0.68 

25.48 

Abell 611 

0.288 

F606W 

-0.39 

-0.11 

0.82 

25.71 

MS 2137-2353 

0.313 

F606W 

-0.27 

0.03 

0.89 

25.91 

RXJ1532+30 

0.345 

F606W 

-0.20 

0.12 

0.97 

26.09 

RXJ2248-4431 

0.348 

F606W 

-0.20 

0.13 

0.97 

26.10 

MACSJ1115+01 

0.352 

F606W 

-0.20 

0.13 

0.98 

26.11 

MACSJ1720+35 

0.391 

F625W 

-0.23 

0.13 

1.07 

26.21 

MACSJ0416-24 

0.396 

F625W 

-0.23 

0.14 

1.09 

26.22 

MACSJ0429-02 

0.399 

F625W 

-0.23 

0.14 

1.09 

26.23 

MACSJ1206-08 

0.440 

F625W 

-0.01 

0.39 

1.19 

26.58 

MACSJ0329-02 

0.450 

F625W 

-0.01 

0.40 

1.21 

26.61 

RXJ1347-1145 

0.451 

F625W 

-0.01 

0.40 

1.21 

26.61 

MACSJ1311-03 

0.494 

F625W 

0.25 

0.69 

1.31 

27.00 

MACS J1149+22 

0.544 

F625W 

0.49 

0.96 

1.41 

27.37 

MACSJ1423+24 

0.545 

F775W 

-0.78 

-0.30 

1.42 

26.11 

MACSJ2129-07 

0.570 

F775W 

-0.73 

-0.24 

1.47 

26.23 

MACSJ0647+70 

0.584 

F775W 

-0.68 

-0.18 

1.50 

26.31 

MACSJ0744+39 

0.686 

F775W 

-0.50 

0.07 

1.70 

26.77 

CLJ1226+3332 

0.890 

F850LP 

-0.89 

-0.20 

2.07 

26.88 



Figure 1. Left: The central regions of Abell 383 (top) and Abell 1423 (bottom). Right: the same regions with all pixels except those 
counted as ICL masked out, equivalent to a surface brightness threshold of 25 mag/arcsec 2 in the rest-frame B-band. Up is North, left 
is East. Images are approximately 2’50” on each side corresponding to 525 kpc for Abell 383 and 575 kpc for Abell 1423. 
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Figure 2. Relation between BCG stellar mass and redshift. 


Figure 4. Cluster mass (M 200 ) vs redshift for CLASH. 
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Figure 3. BCG stellar mass vs cluster total mass (M 200 ), see 
Table |T] for origin of cluster masses. 


cluster galaxies with spectroscopic redshifts lie within 2<r of 
their photometric redshift. Thus we would expect the per¬ 
centage of companions which would be misidentified non¬ 
cluster members (or cluster members which would be erro¬ 
neously excluded) to be 4% of the number of companions 
counted, thus the error on N c is 4% of N c . 


4 RESULTS 

The results for BCG masses, ICL fractions, numbers of BCG 
companions (N c ) and inferred mass growth from their accre¬ 
tion are summarised in Table [5] 


4.1 Results for BCG masses 

We will define the level of correlation for our results using 
the Spearman rank correlation. The Spearman rank coeffi¬ 
cient (SR) is between zero and one with SR=0 for highly 
uncorrelated variables and SR=1 (-1) for highly correlated 
(anti-correlated) variables. The Spearman rank also includes 
a significance of deviation from zero (SDZ) indicating the 
statistical strength of the correlation. SDZ represents the 
probability that a value greater than or equal to the observed 


SR would be calculated under the null hypothesis (of no cor¬ 
relation), and is given by SDZ=SR(^(n — 2)/(1 — SR 2 )), 
where n is the number of degrees of freedom, and SDZ is 
distributed as a Student’s t-distribution with n — 2 degrees 
of freedom. SDZ has values between 0 and 1 where a small 
value for SDZ indicates a high significance for SR. 

Figure [2] shows the relation between the BCG stellar 
masses and redshift, the Spearman rank coefficient for this 
figure is SR= 0.7 and its statistical significance is SDZ= 
10 -3 , indicating a moderately strong correlation with a high 
statistical significance. This figure clearly shows that the two 
highest redshift clusters house the most massive BCGs in 
the CLASH sample. Even without these two high redshift, 
high mass BCGs, the rest of the sample still forms a fairly 
significant correlation (SR=0.6, SDZ=10 -2 ). However when 
the two highest redshift points are excluded the slope of the 
gradient for a simple y = mx + c fit decreases by half, form¬ 
ing a fairly flat relation. Figure [3] clearly shows that there is 
no correlation between cluster mass and BCG stellar mass, 
with a Spearman rank coefficient SR= —0.05 and SDZ= 0.8. 
Given that cluster masses and BCG masses are usually seen 
to be positively correlated (e.g., |Stott et a l. |201 0), this sug¬ 
gests that the almost flat trend seen here is due to a selection 
bias of the CLASH sample, whereby clusters with BCGs of 
fairly narrow mass have been selected. The relation between 
cluster mass and redshift is shown in Figure [f] where no 
correlation is seen (SR=0.1 and SDZ=0.6). The majority of 
the clusters in the sample (18/23) fall over a narrow range 
in cluster mass (6 x 10 14 AfQ — 2 x 10 15 Mq) so perhaps it is 
not surprising to see the BCG masses showing little varia¬ 
tion with redshift (with 18/23 BCGs being in a 5 x 10 n Mg 
mass range). We suggest that the BCG masses recorded for 
CLASH are a result of selection, rather than providing ev¬ 
idence for lack of BCG mass growth or an inverted BCG 
mass growth. The absence of any correlation between clus¬ 
ter mass and redshift means subsequent trends should not 
be affected by the variation of cluster mass with redshift. 


4.2 Results for ICL fractions 

Figure [5] shows the fraction of the total cluster light con¬ 
tained in the ICL for the CLASH sample. As mentioned 
above the higher redshift (z > 0.4) clusters in the sample 
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Figure 5. Fraction of total cluster light contained in the 
ICL below a surface brightness threshold of observed frame 25 
mag/arcsec 2 and above lower limits of 25.5 mag/arcsec 2 (black) 
and 26 mag/arcsec 2 (red); and the fraction of the total cluster 
light contained in the BCG (blue stars). 


have cosmology corrected, observed 25 mag/arcsec 2 surface 
brightness thresholds which are fainter than the 2cr sky noise 
level and are thus excluded from any ICL analysis - we limit 
our ICL analysis to below this redshift. 

A trend of increasing ICL fraction with decreasing red- 
shift is clearly seen, with a strong correlation of SR=—0.99 
and a high significance of SDZ= 10~ 10 (for both black and 
red points). The CLASH sample shows an increase in the 
amount of the total cluster light contained in the ICL of a 
factor of ~ 4 — 5 between z = 0.40 and z — 0.19. 

The 2 <j sky background levels for th e CL ASH data are 
~ 25.5 — 26.4mag/arcsec 2 , (see Section 3.2 1 meaning the 
ICL is only measured above the sky for clusters at redshift 
z < 0.4. As redshift increases, cosmological dimming and 
waveband shift effects cause the observed surface bright¬ 
ness equivalent to the rest frame 25 mag/arcsec 2 to be¬ 
come fainter. The upper limit on surface brightness mea¬ 
sured eventually approaches the level of the sky, making 
the ICL unmeasurable below 25 mag/arcsec 2 at higher red- 
shifts. Even with the limit of z < 0.4 some of the higher 
redshift clusters left over have their 25 mag/arcsec 2 thresh¬ 
old limit within 0.5 mag/arcsec 2 of the sky. In order to test 
the possible bias introduced by measuring the ICL close 
to the sky noise limit we measured the ICL over surface 
brightness ranges of 0.5 and 1.0 mag/arcsec 2 below the 25 
mag/arcsec 2 threshold (25.5 and 26 mag/arcsec 2 ). The re¬ 
sults of this are shown in Figure [5] where the black points 
indicate the range 25-25.5 mag/arcsec 2 and the red points 
indicate the 25-26 mag/arcsec 2 range. We find that a fainter 
lower limit of 26 compared to 25.5 mag/ arcsec 2 increases 
the ICL fraction measured by an average of 41±3% (i.e. in¬ 
creasing the amount of the total cluster light in the ICL from 
~16% to 23%). Lowering the faint end limit above which 
the ICL is measured has a strong effect on the amount of 
ICL recorded. We will discuss the effect of the lower sur¬ 
face brightness limit on the ICL recovered in Section [5] The 
steepness of slope of the ICL fraction with redshift recorded 
for both surface brightness lower limits shows that the ICL 
really is evolving rapidly from z = 0.4 — 0.2. 


Figure 6. The fraction of the total cluster light contained in the 
ICL for CLASH vs BCG stellar mass. 


The higher redshift study of Burke et al. 2012 (hereafter 
B12) examined clusters at 0.8< z <1.2 using near infrared 
data and showed a quadrupling in the cluster’s ICL fraction 
between their average redshift z ~1 and the present. The 
slope of the points in Figure [5] shows a fairly steep increase in 
the ICL fraction between z = 0.40 and z = 0.19, indicating 
a rapid but steady buildup of the ICL between intermediate 
redshifts and the present. This also supports our previous 
conclusions in BC13 that if many mergers occur in the cores 
of clusters, rather than increasing the mass of the BCG, the 
stellar mass in the merging galaxies ends up contributing to 
the ICL, allowing it to grow rapidly at late times. 

Other findings in the literature for the growth of the ICL 
over similar redshifts show wide ranging values for ICL frac¬ 


tions. Zibetti et al. (2005) find SDSS clusters at z = 0.2 —0.3 


to have 11% ICL; Giallongo et al. (2014 1 find 23% ICL in a 


cluster at z = 0.4; Feldmeier et al. (20041 find 8-28% ICL 
below a surface brightness threshold of V=26 mag/arcsec 2 
(similar to that used here) in clusters at z = 0.16 — 0.19 and 
Krick & Bernstein (20071 find 6-22% ICL below the same 


threshold for clusters at z = 0.05 — 0.3 with an increase in 
the ICL fraction from 11% to 14% over the redshift range 
studied. Simulations of ICL growth also show large ICL frac¬ 
tions at z = 0, with typical values in the range of 15-45% 


of the total cluster mass or light in the ICL (Contini et al. 

2014j Purcell et al. 2007; Henriques & Thomas 2010 

Mu- 

rante et al.]|2007 

Puchwein et al. 2010; R.udick et al. 

2011 

Budzynski et al. 

2014), values which are in general good 


agreement with the low redshift end of the CLASH data. 
Both the observational and simulated results from the liter¬ 
ature show no tight consensus on the fraction of total cluster 
light contained in the ICL at low redshifts (0 < z < 0.3), 
with measurements in the general range of 10-30% over a 
range of cluster masses, wavebands and surface brightness 
thresholds. The observed ICL fractions found here fit rea¬ 
sonably well (qualitatively) within the range found in other 
observational studies, even though the clusters and methods 
of measurement vary between studies. However, the lower 
surface brightness limits of the measurements in these other 
studies are not well described, and as discussed above this 
may significantly affect the amount of ICL measured. 

In Figure [6] we compare the BCG masses to the ICL 
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Table 4. Results of partial Spearman rank test for the correlation 
of redshift, BCG stellar mass and ICL fraction. 


Spearman Rank Partial Spearman rank 


Rbcg ICL = -0.7 Rbcg ICL,z = 0.2 
RiCL z = -0.9 RiCL z,BCG= -0.8 

Rbcg z = 0.7 Rbcg z,icl= 0.2 


fractions. This figure shows an increase in ICL fraction with 
decreasing BCG mass, with a fairly strong correlation of 
SR=—0.7, with a high statistical significance SDZ= 0.01. 
However, the clusters with the highest ICL fractions are the 
lowest redshift clusters, and the lowest mass BCGs are also 
at lower redshifts (see Figure [2j, suggesting that the correla¬ 
tion seen is purely due to redshift. A partial Spearman rank 
correlation allows one to determine the correlation between 
two variables A and B, given that they are both also corre¬ 
lated with a third variable C, allowing one to quantify the 
correlation between A and B once their independent corre¬ 
lations with C have been removed. The partial Spearman 
rank coefficient has the form, 


SRab,c = 


SRab — SRacSRbc 


(5) 


[{1-SR\ C ){1~SR% C )]^' 

As with the standard Spearman rank coefficient, a value 
close to ±1 indicates very correlated (or anti-correlated) 
variables and a value close to zero indicates uncorrelated 
variables. For a more detailed discussion of the partial Spear¬ 
man rank see Collins & Mann (1998). Given that the BCG 
mass and ICL fraction are both correlated with redshift we 
carried out a partial Spearman rank analysis for these three 
variables. The results of this analysis are shown in Tabic [4] 
alongside their original Spearman rank coefficient. From this 
table it is clear that the ICL fraction is correlated with red¬ 
shift, independent of the BCG mass, but the ICL fraction 
and BCG mass are not at all correlated with each other 
when the effect of redshift is included. Therefore the signif¬ 
icant correlation found for the ICL fraction and BCG mass 
from the original Spearman rank analysis is actually spuri¬ 
ous. The partial Spearman rank for BCG mass vs redshift 
is now only 0.2 when the strong ICL vs redshift correlation 
is accounted for. 

Figure [7] shows the ICL fraction with redshift when the 
CLASH results are compared to the predicted ICL fractions 
from the simulations of|Rudick et al.| (|2011[). The |Rudick| 


et al. (2011) data points are for ICL measured below a 
surface brightness threshold of 25 mag/arcsec 2 in the V- 
band, which is slightly redder than the observed B-band for 
CLASH. In this case one might expect the simulated clus¬ 
ters to have higher ICL fractions as the ICL is thought to be 
built up from mergers between the red, central galaxies in 
the core, and in past studies the ICL has been shown to have 
an old stellar population (e.g., Krick fe Bernstein|2007 l. 

The slope of the observed and simulated data clearly dif¬ 
fer from each other, with the observed ICL fractions showing 
much steeper increase with redshift than is predicted by the 
simulations. One possible cause for this difference may be 
the depth of the observational data compared to the almost 
unlimited depth of the simulated data. As discussed above 
and illustrated in Figure [5j a change in the lower limit of 
the flux available to be measured can cause a significant 



Redshift 


Figure 7. Comparison of CLASH (black) ICL fraction measure¬ 
ments with the predictions of the simulations of |Rudick et al.| 
(20lTj blue). The dashed line shows the best fit for the simulated 
data and the dotted line shows the best fit for the CLASH data. 


change in the ICL fraction recovered. Since simulated data 
can probe surface brightness depths at faintness levels well 
beyond that of observational data it is possible that simula¬ 
tions will predict significantly more ICL than is observed. In 
the case of Rudick et al. ([201 lj they are able to measure the 
ICL down to surface brightnesses of uv =35 mag/arcsec 2 , 
this level is far below the sensitivity of most current observa¬ 
tional data, and the amount and extent of extra flux which 
may be detected at these faint surface brightness levels is 
currently unknown. 

The simulations of Rudick et al. ( |2011[ ) have cluster 
masses which are generally smaller than those of the CLASH 
sample, approximately by an order of magnitude in most 
cases. However for the order of magnitude range of cluster 
masses studied by Rudick et al.! (2011) there is very lit¬ 
tle variation in the ICL fraction (see their Figure 4) mak¬ 
ing the comparison between these simulations and CLASH 
valid. Other simulations of the assembly of the ICL generally 
show a substantial growth over the redshift range 2 = 1 — 0, 
|Contini et al.|(|2014 1 predict a factor of 6 growth in the ICL 
and Murante et al. (2007) predict factor of 3-5 times growth 
for the ICL, both of which are similar to the growth found 
here (factor of 4-5) and in Burke et al. (|2012, factor of 4 


from 2 = 1.2 to 2 = 0). 

When we compare the fraction of the total cluster light 
that is contained in the ICL with the mass of the host clus¬ 
ters we find a fairly flat distribution, with a large scatter 
in ICL fraction especially at the lower mass end. No strong 
trend is found between ICL fraction and cluster mass (SR= 


-0.1 and SDZ= 0.8). Krick et al. 

(20061 and Krick & Bern- 

stein 

(2007) also find no trend 

oetween cluster mass and 


ICL fraction in their 2 = 0.05 — 0.3 Abell clusters. This lack 
of correlation is also generally found in simulations, where 
there is no dependence between the fraction of total cluster 


light or mass contained in the ICL and the host halo mass 

(Contini et al. 2014 Henriques & Thomas 2010 

Puchwein 

et al. 2010 Rudick et al. 2011), however in the si 

mulations 
a positive 

by Purcell et al. (2007) and Murante et al. (2007) 


correlation between halo mass and ICL fraction is seen. 


The fraction of the total cluster light which is contained 
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Figure 8. Fraction of total cluster light contained in the 
BCG+ICL against cluster total mass. Black points are from 
CLASH, blue crosses are the results of |Gonzalez et al. ([2013]) 
at an average redshift of z = 0.1. 


in the BCG was measured from the same aperture magni¬ 
tudes which were used to determine the stellar mass. The 
fraction of the total cluster light contained within the BCG 
and ICL for the CLASH sample is shown as the blue points 
in Figure [5] From this figure it is clear that the majority of 
the BCG+ICL light is in the ICL for the CLASH clusters 
at redshifts 2 < 0.4, with the ICL containing 70-80% of the 
BCG+ ICL light at these redshifts. This result is in line with 
the findings of previous studies which fit the surface bright¬ 
ness profile of the BCG+ICL (e.g., |Gonzalez et al.||2007[ ). 
At redshift 2 ~ 0.4 the ICL and BCG both contain similar 
fractions of the total cluster light. After 2 < 0.4 the ICL 
grows rapidly but the BCG stays fairly constant in terms 
of its fraction of the total cluster light. This indicates that 
the BCG has assembled the majority of its stellar mass by 
2 ~ 0.4 at the latest, and the major evolution or merging 
events after this point are causing the growth of the ICL 
rather than the BCG. 

The BCG+ICL light as a percentage of the total clus¬ 
ter light is plotted against the total cluster mass in Figure[8j 
with the results of Gonzalez et al. ( 2013| ) for comparison. In 
this figure there is no trend with mass for the CLASH re¬ 
sults on their own (SR= 0.1, SDZ=0.8). Previous studies 
find that lower mass clusters have higher BCG+ICL frac- 


tions (e.g., Gonzalez et al. 

2007, 2013 

Toledo et al.||2011|). 

Over the cluster mass range 0.9 x 10 14 

- 6 x 1O 14 M 0 Gon- 


zalez et al. (2013) report BCG+ICL fractions of 20-50%, 


similar to the BCG+ICL fractions found for the low mass 
end of the present study, however the lower mass points of 


Gonzalez et al. (2013) have a lower average redshift than 


the CLASH clusters of 2 = 0.1 (c.f. 2 = 0.3 for CLASH) 
and therefore may have higher ICL fractions due to having 
had more time for the clusters to grow their ICL. When 
the CLASH and Gonzalez points are considered together 
they form a strong correlation with SR=-0.8, and high sig¬ 
nificance, SDZ= 10~ 6 . This correlation could be due to the 
changing fraction of ICL between clusters of different masses 
or it could be due to the evolution of the ICL with redshift. 



Redshift 



Redshift 


Figure 9. Top: The numbers of BCG companions ( N c ) within 
a 50 kpc aperture centred on the BCG against redshift. Bottom: 
Total mass of BCG companions within 50 kpc as a function of 
redshift. For both plots black points show raw total companion 
numbers or masses, red crosses show just those companions which 
will merge with the BCG according to their dynamical friction 
timescale. 


4.3 Results for BCG companions and mergers 


Now we turn our attention to the numbers of BCG compan¬ 
ions. The top panel of Figure [9] shows the redshift relation 
of the total number of BCG companions within 50 kpc and 
the number of these companions for which their dynami¬ 
cal friction timescale indicates that they will merge with 
the BCG in the time between the cluster redshift and the 
present. Both the total N c and the merging N c are highly 
uncorrelated with redshift, with SR= —0.1 (SDZ= 0.6) and 


SR=0.45 (SDZ= 0.03) respectively. Figure 10 shows the re¬ 


lation between N c and redshift when the number of BCG 
companions in CLASH is limited to the same magnitude 
completeness as that of BC13. This figure shows a decrease 
in the number of BCG companions with redshift, however 
this trend is fairly weak, with SR=0.5, but highly significant, 
all the companions have luminosi- 


10 


SDZ=10 . I 11 Figure _ 

ties greater than 1:20 of the BCG, whereas Figure [9] shows 
companions down to a luminosity ratio of 1:1000 of the BCG 
for the total N c - thus this implies that for more massive BCG 
companions there is a general decrease in their numbers with 
redshift. This decrease is certainly to be expected for the 
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Figure 10. The number of BCG companions (Nc) with redshift 
from CLASH (black) and Burke & Collins 2013 (red crosses), 
where the detection limit for companions is consistent between 
the two samples. The errors on the CLASH points are smaller 
than those on the BC13 points due to the more accurate method 
by which non-cluster galaxies were excluded from the companion 
counts, see Section [3] 


most massive BCG companions as they have much shorter 
merging timescales than less massive companions, and this 
shows that there are several major and large-sized minor 
mergers (BCG:companion masses ~ 1:3-1:10) in the cluster 
core between 2 < 1.5 and the present. This also suggests 
that the larger total N c shown in the top panel of Figure [9] 
is a result of the presence of many low luminosity companion 
galaxies, which are not detected in BC13. The presence of 
more low luminosity galaxies is corroborated in Figures [15] 
and[l6]and is discussed below. It is worth noting that for the 
one overlapping cluster between the two samples, CLJ1226, 
the same number of companions is measured from the data 
used here as from BC13 when the same magnitude limit for 
companion detection is applied for both datasets. 

The number of BCG companions shows no clear trend 
with cluster mass, for CLASH companions that both will 
and won’t merge (top panel Figure [IT] black points SR=0.4, 
SDZ= 0.08; red points SR=0.4, SDZ=0.05); clear trends are 
also not seen when the higher redshift BC13 sample is in¬ 
cluded (SR=—0.2, SDZ=0.2). The top panel of Figure [15] 
shows the relation between the number of BCG companions 
within 50 kpc and BCG stellar mass. A possible weak trend 
is found in this plot (SR=—0.53 —0.53, SDZ= 10 -3 ) but 
no trend is seen for the number of companions which will 
merge (SR=0.22, SDZ=0.3). Given that CLASH represents 
a sample with a small range of BCG masses it is perhaps not 
surprising to see no trend in the numbers of merging BCG 
companions when compared to BCG masses. 

The bottom panels of Figures |9j [TT] and [12] show the 
total stellar mass contained in the BCG companions com¬ 
pared to redshift, cluster total mass and BCG stellar mass 
respectively. The bottom panel of Figure [9] shows both the 
total companion stellar masses and the stellar masses of the 
companions which will merge with the BCG by 2 = 0. This 
figure shows that the majority of the total companion stellar 
mass, on average 68% of the total, will merge with the BCG. 
However, in the top panel of this figure it is shown that only 



Cluster mass / Mo 



Cluster mass / Mo 

Figure 11. Top: The relation between the number of BCG com¬ 
panions within 50 kpc and cluster mass (M2oo)- Black points 
show total companions, red crosses show companions which will 
merge with the BCG. Bottom: Total mass of BCG companions 
which will merge with the BCG according to dynamical friction 
timescale, compared to the cluster total mass. 


around 30% of the total number of companions will merge. 
It seems that the merging companions are indeed more mas¬ 
sive, and the mass contained in all the less massive compan¬ 
ions is smaller despite their larger numbers. However, there 
is no correlation between redshift and total companion stel¬ 
lar mass both for all the companions (SR=0.1, SDZ= 0.7) 
and just those that will merge (SR=0.2, SDZ=0.4). 

The bottom panels of Figures E3 and [l2] show only the 
mass in the companions which will merge with the BCG 
by 2 = 0. Figure [12] shows no correlation between BCG 
mass and merging companion mass (SR=—0.05, SDZ=0.8) 
and Figure [TT] shows a weak, positive correlation between 
total cluster mass and merging companion mass (SR=0.6, 
SDZ= 0.01) suggesting the possibility that more massive 
companions reside in more massive halos. In a recent study 
by |Tal et~ al. ( 2014) of massive galaxies in groups at 0.2< 2 < 
1.2 it was found that in more massive haloes, more of the 
total group mass was in the satellites of the central galaxy. It 
was also found that the companion masses are not correlated 
with redshift for fixed central galaxy mass, consistent with 
what is seen here. Tal et al. (2014) find that more massive 
central galaxies have more massive companions which is not 
seen here, however CLASH shows no correlation between 
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Figure 12. Top: The number of BCG companions (Nc) within 
50 kpc vs BCG mass for CLASH. Black points show total com¬ 
panions, red crosses show companions which will merge with the 
BCG. Bottom: Total mass of BCG companions within 50 kpc 
that will merge with the BCG according to dynamical friction 
timescale, compared to the BCG mass at its redshift. 


BCG mass and halo mass - a trend which is usually observed 
- so the absence of this trend here is most likely due to the 
small dynamic range of BCG masses. 

When we examine the relation between the number of 
companions to the BCG and the fraction of the total cluster 
light contained in the ICL (for both the total N c and just 
those that will merge) we find a large scatter and no obvious 
trends. The total N c has no correlation with ICL fraction 
with SR=0.01 (SDZ=0.9), and the numbers of merging com¬ 
panions also show very little evidence of correlation with the 
ICL fraction (SR=—0.4, SDZ=0.03). 

The merging companion masses also show no trend 
when compared to the ICL fractions for each cluster 
(SR=—0.08, SDZ=0.7). When the BCG mass and compan¬ 
ion masses are summed together there is also very little ev¬ 
idence for correlation ( SR=—0.6, SDZ=0.02). This lack of 
correlation suggests that the numbers of companions and 
their masses has at best a weak effect on the amount of ICL 
at a given redshift. However the number of companions at 
a given redshift does not provide any information about the 
number of mergers that the cluster core may have experi- 


o 



Figure 13. The inferred mass growth additional to the BCG 
mass measured at the cluster redshift from the accretion of its 
companions within a 50 kpc aperture. Top: black points show 
the measured BCG mass at its current redshift, red points show 
the BCG mass after the accretion of its companions. Bottom: the 
mass growth of the BCG from its accreted companions divided 
by the BCG’s original mass. 



Figure 14. BCG mass vs the BCG mass after merging of com¬ 
panions with dynamical friction timescales less than the time be¬ 
tween z = 2 (BCG) and 2 = 0. The solid line represents x=y, the 
dashed line is the best fit for a power-law fit. 


enced, hence this gives no insight into the effect of mergers 
on the growth of the ICL fraction. 

Figure [13] shows the BCG masses after the accretion of 
their companions as a function of redshift, with the orig¬ 
inal BCG mass for comparison. The BCG masses tend to 
be similar before and after the addition of the merged com¬ 
panion masses in most cases. For a linear fit the gradients 
for both sets of points are within 2cr of each other, showing 
no redshift dependance on final BCG mass. The lower panel 
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of Figure [13] shows the BCG mass growth divided by the 
original BCG mass against redshift, and this also shows no 
obvious trend in mass growth with redshift. For 11 of the 
23 clusters the BCG mass after accretion of companions is 
the same within the errors as the BCG mass at the cluster 
redshift. Figure [14] shows this more clearly, where the BCG 
measured mass is plotted against the BCG mass after accre¬ 
tion of its companions. The solid line on this figure indicates 
y = x and the dashed line is a geometric y = ax° + c best fit, 
for which a=(3±2) x 10 -4 , b=1.29±0.2, c=(1.93±0.3)xl0 11 . 
The possible slight increase in the steepness of the curve ap¬ 
pears to be from the growth of the least massive BCGs; 
however both the change of the slope and the mass growth 
of the BCGs are moderate at best. The maximum stellar 
mass growth for the whole sample is 3.4 times the original 
BCG mass, but the mean growth is more moderate at 1.4 
times. The amount of accreted matter from a merger which 
ends up falling on to the central BCG or is instead dis¬ 
tributed within the diffuse ICL is observationally undefined. 
Simulations of the buildup of the ICL predict that the frac¬ 
tion of the stellar matter which goes into ICL rather than 


BCG during a merger is between 30-80% ( 

Murante et al. 

2007 

Laporte et al.|2013| jConroy et al.;|2007; |Contini et al. 

2014 

Puchwein et al. 

2010 

I. Assuming an average of 50% 


of merging stellar mass goes into ICL rather than BCG, the 
subsequent BCG mass growth for the CLASH sample should 
be 1.2 times the original BCG mass. This suggests that on 
average BCGs do not grow by a large amount in stellar mass 
from the accretion of their companions. 

The BCG after accretion of its companions shows no 
correlation with cluster mass (SR=0.2, SDZ= 0.4). The 
BCG growth also shows no significant correlation with red¬ 


shift or original BCG mass (as shown in Figures 13 and 141 
and we find only moderate growth in stellar mass. However, 
the ICL shows no significant increase over the same redshift 
range and it is apparent that the growth of the ICL is the 
major evolution event in the cluster core during the second 
half of the lifetime of the Universe. 

Figure [15] shows the contribution to the stellar mass 
growth of the BCGs that accreted companions will make as 
a function of their stellar mass ratio with respect to the BCG 
mass. The CLASH results are shown for companions which 
will be accreted (black) and for companions as limited to 
the same completeness as BC13 (red), alongside the results 
for BC13 (blue). This figure shows that ~25% of the total 
merging mass is in companions of mass ratio to the BCG of 
1:5-1:10 and 1:10-1:20. A further ~16% comes from com¬ 
panions of mass ratios of both 1:2-1:3 and 1:3—1:5 with the 
BCG. Finally only ~6% of the merging companion mass is 
in galaxies of 1:20-1:100 of the mass of the BCG. Added 
together the majority of the merging mass in BCG compan¬ 
ions, around 60% is in low mass BCG companions of masses 
less than 1/5 of the BCG, showing that for CLASH the ma¬ 
jority of mergers in the cluster core will be minor mergers 
and there will be very few major mergers. When compared 
with the results of BC13 (for the points which are limited 
to the same companion brightness completeness limit) it is 
clear that there is an abundance of low mass companions in 
the CLASH cluster cores compared to the sample used in 
BC13. This is to be expected as the BC13 sample is at a 
higher average redshift than CLASH (c.f. z = 1 for BC13 
and z = 0.4 for CLASH) and the more massive companions 



Mass ratio (BCG:companion) 


Figure 15. The fraction of the BCG mass growth from the ac¬ 
cretion of its companions between the cluster redshift and the 
present compared to the mass of the companion - companion mass 
decreases from left to right. Black: total merging companions in 
CLASH; red: CLASH companions limited to the same luminosity 
as the completeness of Burke fe ColIins| [2013); blue: results of 
|Burke fc Collins] (|2013[) ■ 



Luminosity ratio (BCG : companion) 


Figure 16. The numbers of BCG companions split by luminosity 
ratio to the BCG compared to results from the literature. Black: 
this study, red: BC1 3, blue: |Edwards Sz Patton] ( |2012| z ~ 0.3), 
green: simulations of|Laporte et al.|(|2013| z =2—0). 


which have shorter dynamical friction timescales will have 
merged in the cluster core between the average sample red- 
shifts of a = 1 and z = 0.4. This indicates that BCGs at 
higher redshifts are built by major mergers (or that major 
mergers are the dominant source of stellar mass assembly in 
the cluster core), and minor mergers are the more dominant 
merging path at lower redshifts. 

When the companions are split into luminosity ratio 
(and hence mass ratio) bins, as would be expected there 
are fewer BCG companions with masses comparable to the 
BCG, and there are very many small mass companions at all 
redshifts. Figure [16] shows the mean numbers of BCG com¬ 
panions by luminosity ratio bin compared to results from 
previous studies. The CLASH clusters show numbers of com- 
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panions comparable to (the same within the errors) those 
found by Edwards fc Patton] (2012), who examined a sample 
of clusters at z ~0.3. The mean redshift of the CLASH sam¬ 
ple is 2 = 0.4, so it is not surprising to find good agreement 
between these results. The other two studies included in this 


figure are the observational study of Burke & Collins 


and the simulated predictions of Laporte et al. (2013), both 


(2013) 


of which examine BCGs at higher redshifts, 0.8 < z < 1.2 
and z ~2 respectively. Since less time has passed for the 
higher redshift clusters it is not surprising that these higher 
redshift studies show and predict a larger number of com¬ 
panions in all mass ratio bins, and in agreement with the 
results shown in Figure [10] there is an overall trend between 
redshift and the total number of BCG companions between 
2 ~ 1.5 and z = 0, which is more obvious when the BCG 
companions are split by their comparative luminosity ratio 
to the BCG. 


5 DISCUSSION 

Across the CLASH sample the stellar masses of the BCGs 
and the total cluster masses both cover a fairly narrow mass 
range. This is also reflected in the lack of correlation between 
BCG stellar mass and cluster mass, where one would nor¬ 
mally expect to see a positive correlation between the two. 
This narrow mass range of both clusters and BCG stars 
is likely to be a selection effect. However, given this fairly 
consistent mass sample, we find strong evidence for the sub¬ 
stantial growth of the ICL and evidence of mergers which 
would lead to a moderate growth of BCG stellar mass or 
alternatively an increase in stars in the ICL. 

Typically galaxy clusters which exhibit strong lensing 
features are extended along the line of sight. If this were the 
case we would expect to measure higher than ‘normal’ ICL 
fractions and numbers of BCG companions for the strong 
lensing CLASH clusters. The selection criteria as described 
in 


Postman et al. (2012) was specifically designed to avoid 


this, with clusters selected for their spherical symmetry in X- 
rays. However if these clusters are highly relaxed as their X- 
ray symmetry would suggest then it may be that we should 
expect them to be highly evolved, having already undergone 
the majority of their galaxy merging interactions thus show¬ 
ing high BCG masses, large ICL fractions and low numbers 
of BCG companions or only low-mass BCG companions. 
Eight of the clusters in the sample are reported to show some 
small degree of substructure (A209, MACS0329, CLJ1226, 
MACS0744, MACS1206, RXJ1347, A2261, RXJ2248), mak¬ 
ing them not fully relaxed. These eight are spread across the 
whole redshift range of CLASH and do not show any notice¬ 
able deviation from the general trends seen in the results, 
suggesting that selection criteria is not a major issue for the 
results from this study. 

According to the dynamical friction timescale for merg¬ 
ers of BCG companions, the BCG stellar masses only in¬ 
crease moderately over the time between 0.18^ z ^0.9, on 
average growing by a factor of 1.4. The eventual resting place 
of accreted matter in cluster cores is unknown, and when 
we assume 50% of merging mass ends up in ICL rather than 
centrally on the BCG, in line with the predictions of sim¬ 
ulations, the mass growth of the BCGs is only a factor of 
1.2 on average. Given the observed large growth of the ICL 


over the same time, it seems likely that the overwhelming 
majority of companion mass must end up in ICL rather than 
BCG. 

The growth factor for BCGs from mergers of 1.4 (1.2 if 
we assume 50% of accreted matter ends up in the ICL) re¬ 
ported here is fairly consistent with the stellar mass growth 
estimate of 1.8 for BCGs in the time between z = 0.9-0.2 


reported by Lidman et al. (20121. Collins et al. (20091 find 
>90% of BCG stellar masses in place by 2 = 1; and Stott | 


et al. (20101 find BCGs at 2 = 1 on average have 95% of the 


stellar mass of those measured locally, which corresponds to 
growth of factors of 1.1 and > 1.05 times respectively. This is 
also fairly consistent with the factor of 1.4 or 1.2 found here. 


Laporte et al. (2013) simulate BCG mass growth by mergers 


and predict BCG stellar mass growth factors of 1.5 between 
0 < 2 < 0.3, 1.9 between 0.3 < 2 < 1 and 2.6 between 
0 < 2 < 1. They also predict that ~30% of accreting mass 
goes into the ICL rather than the BCG. The simulation of 


De Lucia & Blaizot (20071 predict BCG stellar mass growth 
of a factor of 4 since 2 ~1 and a factor of 1.6 since 2 ~0.4, 


however De Lucia & Blaizot (2007) have no prescription for 


the ICL in their simulation. The implied mass growth we 
find for CLASH lies in between the previous observational 
results of Collins et al., Stott et al. and Lidman et al., but 
is notably less than the growth predicted by simulations. 

The stellar population of the ICL is largely unknown, 
and if it is built from the accretion of other cluster galaxies 
then its stellar population may be quite varied from loca¬ 
tion to location within a cluster. Available estimates use the 


colour of the ICL to derive its stellar population (e.g., Krick 
& Bernstein 2007 \ and suggest that the ICL near the clus¬ 


ter core has a population similar to that of BCGs, being 
old and red (however Krick fe Bernstein|2007 also find gra¬ 
dients in their ICL colours). If we assume that the stellar 
population of the ICL is the same as that of the BCG we 
can give a rough estimate of its mass and its mass growth 
over the redshift range of this study. Figure [5] shows the rel¬ 
ative contributions of the BCG and ICL to the total cluster 
light. If we assume that the stellar mass scales directly with 
stellar luminosity for the same kind of stellar population, 
this figure would suggest that the ICL has 5-10 times more 
stellar mass than the BCG at low redshifts 2 < 0.4, and a 
similar mass to the BCG at higher redshifts 2 ~ 0.4. The 
average BCG mass for CLASH is ~ 4 x 10 n Mg. Taking 
a conservative value of 5 times the mass of the BCG gives 
an average estimate for the ICL stellar mass at 2 < 0.4 of 
2x10 12 Mq. Given this estimate for the ICL mass at low 
redshift, the ICL must accrete approximately 1.6x10 12 Mq 
of stellar mass, presumably from merging companions since 
2 ~ 0.4. This is a significantly higher mass than is contained 
in all the companions to the BCG within 50 kpc for all the 
CLASH clusters, by around a factor of 5-10. If the mass in 
the ICL is indeed this high, and it grows substantially via 
tidal stripping from mergers, then the majority of the mass 
which grows the ICL must come from outside of 50 kpc, as is 


predicted by some simulations of ICL growth (Conroy et al. 
2007[ [Murante et al.|[2007[ |Purcell et al.||2007| |Puchwein 


et al. [2010). 


An amount of stellar mass accretion of the magnitude 
estimated above adds up to a lot of mergers, especially if 
the amount of accreted matter that ends up deposited on 
the ICL is only 50% of the total merging mass. This sug- 
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gests that our assumed value of 50% may be a very con¬ 
servative estimate, and that the actual fraction of accreted 
matter from mergers which contributes to the ICL may be 
significantly higher. 

We find that the fraction of the total cluster light which 
is contained within the ICL grows significantly, by a factor 
of 4-5 since z = 0.4. This result is broadly in line with pre¬ 
vious observational studies and is consistent with the pre¬ 
dictions of simulations which make special consideration for 
the ICL, however the rate at which the ICL grows is some¬ 
what different to that predicted by simulations. Given that 
simulations generally predict the ICL to assemble through 
galaxy mergers and interactions with the BCG (including 
tidal stripping), it is clear that the consideration of the ICL 
is vital when examining the growth and assembly of BCGs, 
especially in light of the, at most, moderate mass growth 
observed for BCGs. 

We find that the amount of ICL detected is not only de¬ 
pendant on the upper limit of the surface brightness thresh¬ 
old below which it is measured, but it is also strongly depen¬ 
dant on the faint-end surface brightness limit. As is shown 
in Figure [5] we find a 40% difference in the fraction of ICL 
recovered when we change the faint surface brightness limit 
by 0.5 mag/arcsec 2 . The faint-end, lower surface brightness 
limit above which flux is measured is generally not discussed 
in previous studies of ICL, both observational and simulated. 
To illustrate this, Presotto et al. (20141 examine the ICL 
in one of the clusters in CLASH, MACSJ1206 (z = 0.44), 
using data from the Subaru Suprime-Cam. They measure 
the ICL to contain 12.5±0.6% of the total cluster light be¬ 
low a surface brightness threshold of fj,y =26.5 mag/arcsec 2 . 
The Presotto et al. (2014) data is deeper than the CLASH 
HST data used here and therefore can measure ICL to a 
rest-frame surface brightness nv =28.5 mag/arcsec 2 , signif¬ 
icantly deeper than is possible here. In order to compare 
more accurately between different observational studies and 
between observations and simulations of the ICL the effect 
of the surface brightness lower limit needs to be more fully 
quantified. 

When we examine the respective contributions of the 
BCG and ICL to the total cluster light we find that the 
ICL contains 70-80% of BCG+ICL light at redshifts a < 
0.2 — 0.3. The fraction of BCG+ICL light in the ICL shows 
a significant increase as redshift decreases. This result also 
points to the importance of the ICL as a component within 
cluster cores and again shows the significant growth of the 
ICL in contrast to the only moderate growth of the BCG. 
The ICL having more of the total cluster light than the BCG 


is also observed in previous studies (e.g. Gonzalez et al. 


2007 20131, however these studies use lower cluster-mass 


and redshift samples than CLASH. The CLASH clusters 
on their own do not show the negative correlation between 
BCG+ICL light as a fraction of the total cluster light and 
cluster total mass as is observed in these previous studies 
(see Figure [8]), however this may be due to the smaller dy¬ 
namic range of cluster masses in CLASH or the lower red- 
shifts of the samples used in other studies. 

When comparing the numbers of companions between 
studies across a range of different redshifts we see fairly 
strong evidence for the occurrence of mergers in cluster 
cores. The BCG companions within a 50 kpc radius are gen¬ 
erally less massive at lower redshift and there are generally 


fewer of them. Similar to the comparable redshift study of 


Edwards & Patton (2012 I, we find a greater number of small- 


mass companions (mass ratios 20:1) than large-mass com¬ 
panions, and find very large companions (masses ^2:1 w.r.t. 
BCG) to be rare. However at the higher redshifts probed by 


Burke & Collins (20131, significantly more large mass com¬ 


panions are seen. This suggests that major mergers are an 
important mechanism for stellar mass build up at high red¬ 
shift and minor mergers are more dominant at low redshift. 
Given the lack of BCG stellar mass growth coincident with 
a large ICL growth it seems that minor mergers build up 
ICL at medium to low redshift. 

|Edwards fc Patton| ( |2012| ) estimate a factor of 1.1 
growth in BCG stellar mass from the accretion of their 
measured companions since their average redshift a = 0.3 
(10% increase in mass), which is similar to the value of 1.4 
(1.2) found here, however Edwards & Patton (20121 measure 
much lower numbers of companions compared to that found 
for the CLASH sample, with their sample having 0.7 com¬ 
panions per BCG on average, c.f. CLASH has 12.7 compan¬ 
ions per BCG on average, of which 2.4 will merge according 
to dynamical friction timescales. Other observational studies 
of low redshift BCGs also find small numbers of companions 
and major merging to be rare. For example Liu et al. (2009) 
use a sample 515 BCGs at z < 0.2 and find large compan¬ 
ions (with luminosity ratios less than 1:4) within 30 kpc of 
the BCG to be rare, with an average 0.1 of these compan¬ 
ions per BCG (c.f. 0.7 for CLASH). Liu et al. find 18/515 
of their BCGs to be undergoing major mergers and suggest 
that BCGs should have increased their masses by a factor of 
1.15 (15%) from major dry mergers since z = 0.7. |McIntosl i 
et al.| ( f2008| ) find only 38/845 BCGs at z^0.12 to be under¬ 


going major mergers (0.04 mergers per BCG) and predict 
that massive haloes are growing by 1 -9% per Gyr by major 
mergers (factor Sjl.09 per Gyr). Whilst these studies predict 
similar mass growths to that found for the CLASH sample, 
we do not see any major mergers currently taking place in 
CLASH (this is likely a result of the selection of the sam¬ 
ple for relaxed-looking clusters) and find that minor mergers 
are far more likely at these low redshifts. The higher-redsliift 
sample presented in BC13 generally finds more larger BCG 
companions, suggesting that major merging is more impor¬ 
tant for the stellar mass build-up of high-redshift BCGs than 
nearby BCGs. For companions more massive than 1:5 of the 
BCG, BC13 find 1.3 companions per BCG, twice as many 
as the 0.7 found for CLASH, and as is evident in Figure p~6| 
BC13 find more than twice as many companions in all mass 
bins compared to CLASH, showing that the majority of 
companions of all masses should indeed have merged with 
the BCG in the time between « ~1 and s ~ 0.4. 


6 CONCLUSIONS 

We have studied the BCG stellar mass growth, ICL growth 
and numbers of mergers onto BCGs for a sample of mas¬ 
sive clusters across a redshift range spanning the second half 
of the lifetime of the LTniverse. We have found the clusters 
in the sample to have a narrow range of cluster and BCG 
masses, and that the BCG masses should only grow by a 
small amount from mergers of the surrounding companions, 
on average a factor of 1.4, or a factor of 1.2 when the amount 
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of merging mass which simulations predict will end up in the 
ICL is considered. Conversely we find that the ICL grows 
by a substantial amount over the same timespan and for the 
same clusters, increasing its contribution to the total clus¬ 
ter light by a factor of 4-5 over this time. We also suggest 
that the growth of stellar mass in the ICL is larger than can 
be provided be the BCG companions within 50 kpc, and 
that the majority of the ICL mass must come from galaxies 
which fall from outside of the core of the cluster. One caveat 
to quantifying the ICL growth is the depth at which observa¬ 
tional data is taken and consistency in comparison between 
ICL surface brightness limits used in model predictions. 

We conclude that BCGs must have assembled the ma¬ 
jority of their stellar masses very early on in the life of the 
Universe, before a redshift of z = 1. We suggest that the 
majority of stellar mass growth of BCGs at high redshifts 
(z > 1) is by major merging, and the majority of merging 
activity in the latter half of the lifetime of the Universe is by 
stripping and minor merging in which the majority of the 
accreted or disrupted stellar mass ends up contributing to 
the ICL rather than the BCG. 
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Table 5. Results for CLASH sample BCG stellar masses, fraction of the total cluster light contained in the ICL, total number of companions to BCG within 50 kpc and the number of 
these which will merge with the BCG according to their dynamical friction timescale, the masses of the merging companions and the merging companion mass compered to the BCG 
mass. 


Cluster 

Redshift 

Cluster mass M 200 
(x1O 15 M 0 ) 

BCG mass 
(xlO n M 0 ) 

% of total cluster light 
ICL BCG 

N c total 

N c merge 

Merging companion 
mass (xlO 11 Mq) 

Total companion mass 
/BCG mass 

Abell 383 

0.187 

1.04+0.07 

1.51+0.24 

23.05+0.73 

2.89+1.19 

20.0+0.8 

3.0+0.2 

1.98+0.16 

1.31+0.13 

Abell 209 

0.206 

1.17+0.07 

2.0+0.16 

17.11+0.77 

2.86+1.37 

16.0+0.6 

0.0+0.2 

0+0 

0 

Abell 1423 

0.213 

1.20+0.59 

1.96+0.11 

17.97+0.73 

2.96+1.39 

9.0+0.4 

2.0+0.2 

1.33+0.11 

0.68+0.05 

Abell 2261 

0.224 

1.76+0.18 

1.74+0.18 

16.64+0.78 

1.78+1.16 

11.0+0.4 

4.0+0.3 

4.14+0.33 

2.38+0.15 

RXJ2129+0005 

0.234 

0.73+0.18 

2.21+0.32 

12.89+1.11 

3.04+1.43 

16.0+0.6 

1.0+0.1 

1.16+0.09 

0.52+0.09 

Abell 611 

0.288 

1.03+0.07 

3.45+0.42 

13.02+0.23 

3.49+1.05 

8.0+0.3 

1.0+0.1 

0.31+0.03 

0.09+0.06 

MS 2137-2353 

0.313 

1.26+0.06 

8.38+0.42 

7.58+1.65 

1.05+0.98 

8.0+0.3 

0.0+0.2 

0+0 

0 

RXJ1532+30 

0.345 

0.64+0.09 

2.20+0.05 

7.44+0.17 

2.68+1.13 

13.0+0.5 

1.0+0.1 

0.32+0.03 

0.15+0.02 

RXJ2248-4431 

0.348 

1.40+0.12 

5.33+0.61 

6.44+0.10 

2.42+0.63 

12.0+0.5 

3.0+0.2 

1.31+0.10 

0.24+0.07 

MACSJ1115+01 

0.352 

1.13+0.10 

2.19+0.17 

6.11+0.29 

1.39+1.09 

14.0+0.6 

1.0+0.1 

0.51+0.04 

0.23+0.05 

MACSJ1720+35 

0.391 

0.88+0.08 

3.83+0.94 

3.16+0.06 

1.49+0.31 

11.0+0.4 

2.0+0.2 

0.78+0.06 

0.20+0.13 

MACSJ0416-24 

0.396 

2.50+0.50 

2.89+0.99 

2.69+0.10 

1.08+0.10 

15.0+0.6 

2.0+0.2 

1.50+0.12 

0.52+0.19 

MACSJ0429-02 

0.399 

0.96+0.14 

4.68+0.23 

2.36+0.13 

1.21+0.59 

8.0+0.3 

4.0+0.3 

1.61+0.13 

0.35+0.04 

MACSJ1206-08 

0.440 

1.00+0.11 

4.93+1.39 


1.21+0.12 

12.0+0.5 

3.0+0.2 

1.10+0.09 

0.22+0.15 

MACSJ0329-02 

0.450 

0.86+0.11 

3.69+0.21 


1.73+0.19 

15.0+0.6 

2.0+0.2 

0.79+0.06 

0.21+0.04 

RXJ1347-1145 

0.451 

1.35+0.19 

3.68+0.85 


0.71+0.10 

14.0+0.6 

3.0+0.2 

1.30+0.10 

0.35+0.13 

MACSJ1311-03 

0.494 

0.53+0.04 

4.49+1.05 


1.00+0.13 

9.0+0.4 

1.0+0.1 

0.23+0.02 

0.05+0.12 

MACSJ1149+22 

0.544 

5.10+1.90 

3.05+0.58 


0.53+0.10 

19.0+0.8 

2.0+0.2 

1.10+0.09 

0.36+0.11 

MACSJ1423+24 

0.545 

0.65+0.11 

4.95+0.13 


2.34+0.45 

11.0+0.4 

3.0+0.2 

0.95+0.08 

0.19+0.02 

MACSJ2129-07 

0.570 

3.50+3.10 

2.44+0.17 


0.69+0.21 

18.0+0.7 

4.0+0.3 

1.74+0.14 

0.72+0.06 

MACSJ0647+70 

0.584 

6.80+1.40 

5.13+1.79 


1.00+0.26 

11.0+0.4 

4.0+0.3 

10.43+0.83 

2.03+0.26 

MACSJ0744+39 

0.686 

0.79+0.04 

8.20+0.43 


1.18+0.16 

5.0+0.2 

2.0+0.2 

0.65+0.05 

0.08+0.03 

CLJ1226+3332 

0.890 

1.72+0.11 

15.00+1.99 


1.47+0.09 

15.0+0.6 

8.0+0.6 

6.71+0.54 

0.45+0.08 



